The Riesz fractional advection-diffusion is a result of the mechanics of chaotic dynamics. It’, s of preponderant importance to solve this equation numerically. Moreover, the utilization of Chebyshev polynomials as a base in several mathematical equations shows the exponential rate of convergence. To this approach, we transform the interval of state space into the interval [−, 1, 1] × [−, 1, 1]. Then, we use the operational matrix to discretize fractional operators. Applying the resulting discretization, we obtain a linear system of equations, which leads to the numerical Solution. Examples show the effectiveness of the method.